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ABSTRACT 

Context. The Ca ii H line is one of the strongest lines in the solar spectrum and provides continuous information on the solar atmo- 
sphere from the photosphere to the lower chromosphere. 

Aims. We describe an inversion approach that reproduces observed Can H spectra assuming local thermal equilibrium (LTE). 
Methods. We developed an inversion strategy based on the SIR code that reproduces Ca ii H spectra in the LTE approximation. The 
approach uses a two-step procedure with an archive of pre-calculated spectra to fit the line core and a subsequent iterative modification 
to improve the fit mainly in the line wing. Simultaneous spectra in the 630 nm range can optionally be used to fix the continuum tem- 
perature. The method retrieves one-dimensional (ID) temperature stratifications neglecting lateral radiative transport. Line-of-sight 
velocities are included post-facto by an empirical approach. 

Results. An archive of about 300.000 pre-calculated spectra is more than sufficient to reproduce the line core of observed Can H 
spectra both in quiet Sun and in active regions. The subsequent iterative adjustment of the thermodynamical stratification matches 
observed and best-fit spectra to a level of about 0.5 % of in the line wing and about 1 % of in the line core. 
Conclusions. The successful application of the LTE inversion strategy suggests that inversion schemes based on pre-calculated spectra 
allow one a reliable and relatively fast retrieval of solar properties from observed chromospheric spectra. The approach can be easily 
extended to an ID non-LTE (NLTE) case by a simple exchange of the pre-calculated archive spectra. Using synthetic NLTE spectra 
from numerical three-dimensional (3D) simulations instead will finally allow one to extend the approach from the static ID-case to 
dynamical atmosphere models including the complete 3D radiative transport. 
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1. Introduction 

"The properties of the solar chromospheric plasma are so poorly 
known that a ll careful stu dies would be enlightening, no matter 
how cruder (ICramlll977h 

The spectral lines of Can H at 396.8 nm and Can K at 
393.3 nm are among the strongest and deepest spectral lines in 
the visible solar spectrum. The intensity radiated away inside 
these lines comes from layers between the continuum forming 
layer and the lower chromosphere, spanning a rang e of about 
1-2 Mm of height in the s olar atmosphere ( Vernazza et al. 
119811: ICarlsson & SteinlfT997h . Both lines have thus been used 
intensively for studies of the chromosphere, where useful 
spectral lines are scarce. Ground-based observations are mainly 
limited to Can H and K (e.g., Jensen & Orrall 1963; Beckers 
'1968"; "Liu '1974; 'Teplitskaia & Firstova 1976; Cram & Dami 
1983: Rutten & Uitenbroek.1991: Lites et al..l999: Rezaei et al 
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iTeplitskaya etall l2QQ9i), Ha (e.g. , 



Brav '1973; "Schmieder et al." "1984"; "Tsiropoula et al." "1993"; 
Balasubramaniam et al. 2004; Lopez Ariste et al. 20^5; 
Sanchez- Andrade Nuno et al.l 120081: ' ICauzzi et al.l l2009l: 
Bostancil '2011'), the Can IR triplet (e.g.. iLopez Ariste et al 



^bostancil 'zu IT), tne L:aii IK triplet (e.g.. li^opez Ariste et ai 
20011: IS ocas-Navarro 20 051: iPietarila et al.ll2007bl: ICauzzi eTal 



2008al: Vecchio etali i2009l) and Hei at 1083 nm (iLitesI 
19861 : Riiedi et aP 1 19951: ISocas-Navarro & Elmord 



Sanchez- Andrad eNuno et al.] l2007l: iKuckein et al.l 



iFelipe et al...201(jh by the absorption in the Earth's atmosphere. 



whereas from space also Mg ii h and k. La, and other lines in 
the extreme UV are accessible (e.g., Kneer et al. 1981; Bonnej 
19811: iFontenla et aD Il988l: ICarlsson et al.l 1 19971: iCurdt et all 



2010lf 



Despite the fact that the Ca ii H and K lines provide a tremen- 
dous amount of information on all of the lower solar atmo- 
sphere, it turned out to be rath er difficult to extract th is infor- 
mation from observed spectra (iLinsky & Avrett|[T970l) . Part of 
the spectral lines forms in an atmospheric regime where the 
gas density is so low that the assumption of instantaneous lo- 
cal thermal equilibrium (LTE) is no longer valid, because colli- 
sions are so scarce that the energy is not necessarily distributed 
evenly over all possible degrees of freedom. A consistent treat- 
ment of line formation in non-LTE (NLTE) is, however, at the 
limits of the theoretical des cription and the commonly available 



computational resourc es ( Solanki et al.lll99ll: ICarlsson & SteinI 



120021: iLeenaarts et aDr2009: Martmez-Svkora et al.ll2012h . After 
a few attei npts for a direct analysis of Can spectra in quiet 
Sun (QS) (Liu & Skumanich '1974") or active regions (e.g., 
Teplitskaia & Efendieva 1975) in the 70's and 80's, the topic 
has not received much attention in the recent years. Apart 
from the inversion code f or Can K spect ra described in 
iRouppe van der VoortI (I2002h that is based on HE ine & LinskyI 
(Il974ah and the manual inversion method for the line wings of 



^ The references here and before are only intended to cover roughly 
every decade and research group/instrument. 
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Can H and K used in ISheminova et alJ (l2005h and ISheminoval 
('2012'), only the group led by Prof. Teplitskaya seems to have 
continued with the developme nt of analy s is and/or in v ersion 
tools for Can H and K (Gri gorveva et al.l ll99lL l2QQQl I2Q09I: 
iTeplitskaya & Grigoryeva>2Q09l) / 

Many other studies dealt with Can H and K only in the 
forward modeling direction, i.e., creating spectra from a model 
atmosphere to match some average characteristic profile (e.g. , 
Linskv & Av rett 1970|; lShine& Linskv 1974b; Suemoto 19771 
Kneer&Mattia. 19781 iLites & Skumanich.1982: .Solanki et al.l 



but not with the intention to analyze individual spectra 



in a two-dimensional (2D) field of view (FOV). The formation 
of Ca n H spectra was addresse d in numer i cal NLTE simulations 
(e.g., [Rammacher & Ulmschn eideii 1 1 9921: iRammacher & Cunt3 
l2QQ5l:ICarlsson & Steiniil997i) . but again in the forward model- 
ing approach, i.e., the synthesis of Can spectra from the output 
of a numerical model (see also Uitenbroek 201 1). While this pro- 
vides consistent NLTE spectra, it does not allow one to analyze 
observed spectra directly. 

For photospheric spectra, several inversion codes exist that 
allow one to retrieve physical parameters of the solar atmo- 
sphere from observed spectra on an optical depth scale (e.g . , 
iRuiz Cobo & del Toro Iniestalll992HSocas-Navarro et al.ll2QQll) . 
With (a few) additional assumptions, the inversion results can 
then also be converted to an absolute geometrical height scale 
(e.g., [San chez Almeida & Lites 2000; Puschmann et al. 2005, 
l2010al:lBeckll201lh . The geometrical height scale is essential 
for determining many quantities derived from the magnetic field 
topology su ch as the el ectric current d ensity or the helicity 
in sunspots (Puschmann et al. 2010b; Ruiz Cobo & Pu schmann 
i2012) . For observations of chromospheric spectral lines, similar 
inversion codes a re only partially available (ISocas-Navarro et aP 
119981: iTziotzioul 120071) . For a possible application to Can H 
spectra, the only ne w development is the NLTE inversion code 
NICOL E based on [Socas-Navarro et aP (l2000l) that was used 
in, e.g., Socas-Navarro et al ] (I2006h . 'Pietar ila et al.1 (l2007a) and 



Ide la Cruz Rodriguez et al.l (l2Q12h 7 However, it still needs to be 
tested if it can be successfully applied to Ca n H spectra. 

For the analysis of Ca n H spectra one therefore has to rely on 
methods that are (hopefully) insensitive to NLTE eff'ects, using 
observed quantities as direct as possible to derive solar atmo- 
spheric properties (e.g., [Beck & Rammache rl |2Q 1 Ol) . In the first 
two papers of this series (Beck et al. 2009, 2012, A&A, in press; 
BE09 and BE 12 in the following), we investigated the statistics 
of velocity and intensity oscillations in Ca n spectra, using only 
Ca n H for the first study, and also Ca n IR at 854 nm for the latter 
thanks to new observations covering both lines simultaneously. 
The conversion from observed intensities or velocities to energy 
in these two studies necessarily had to rely on some assumptions 
such as the LTE condition that have an impact on the final re- 
sult. A consistent NLTE analysis would be strongly preferred, 
but lacking a suitable method for the automatic treatment of the 
observations with several 10.000 spectra each, we now first de- 
veloped and applied an LTE inversion strategy for Can spec tra 
based on the SIR code (iRuiz Cobo & del Toro Iniestalll992h as 
preliminary step to a possible full NLTE treatment in the future. 
The LTE assumption does not fail instantly in the solar atmo- 
sphere, but results derived assuming LTE will deviate more and 
more from the "reality" the closer one gets towards the very line 
core of Ca n H that forms highest above the continuum. The LTE 
inversion can also be helpful for an NLTE analysis: it both out- 
lines a suitable approach for the inversion procedure and pro- 
vides a reasonable initial model for the NLTE fit that presumably 
will speed up its convergence. 



The observations used are outlined briefly in Sect.[2l whereas 
Sect. [3] describes the stray-light correction. The LTE inversion 
approach is described in Sect. IH Section \5\ shows the result of 
the application of the code to the spectra. The results are dis- 
cussed in Sect.[6l Section [7] provides our conclusions. Appendix 
IaI explains the animation of the inversion results that accompa- 
nies this paper. 



2. Observations 

For the present study, we used the Can H spectra of the QS 
observations on disc centre labeled No. 1 and 2 in BE 12. The 
first observation is a large-area scan taken on 24 July 2006 
that was used in BE09 before and is described in detail there. 
It consists of simultaneous spectroscopy of Can H and spec- 
tropolarimetry at 630 nm obtained with the POlarimetric LIttrow 
Spectrograph (POLIS, Beck et al. 2005) at the German Vacuum 
Tower Telescope (VTT, Schrote r et al.lll9 85) in Izana, Tenerife, 
Spain. The integration time per scan step was 6.6 sees. The im- 
age of the solar surface was scanned in 150 steps of Of! 5 step 
width, covering a total FOV of about 75'' x 70''. The slit width 
was (y.'5, and the spatial sampling along the slit was (T.'S. The sec- 
ond observation is a time- series taken with POLIS on the same 
day that is described in detail in lBeck et all (l2008 b. BE08). The 
integration time was 3.3 sees per scan step and the cadence of 
repeated co-spatial spectra is about 21 sees. 

For comparison with the results for these QS data on disc 
centre we also applied the inversion code to an observation of 
an active region (AR) off' the disc centre at a heliocentric angle 
of about 50°. These data were taken on 8 December 2007 with 
a combinatio n of POL IS and the Tenerife Infrared Polarimeter 
(TIP, Martine z PiUet et al. 1999; Collados et al. 2007). The data 
consist of three scans of 200 steps with a step width of Off 5, 
covering a great part of the active region NOAA 10978. The 
integration time was 3.3 sees per scan step. The rest of in- 
strumental characteristics for POLIS was as above. Overview 
maps o f these observations can be found in lBeck & Rammached 
(120101), their Fig. 16. One of the AR maps is described in de- 
tail in iBethge et al.l (l2012h . All data were taken with real-time 
seeing correctio n by the Kiepenheuer-Ins titute adaptive optics 
system (^KAOS, von d er Liihe et al.l '2003). The spatial resolu- 
tion of the observations, being slit- spectrograph data, cannot be 
improved post-facto by im age reconstruction techniques (e.g., 
'Puschmann & Beckl l201 ll) apart from a possible correction for 
the instrumental point spread function (.Beck et al...201L BE 11 
in the following). 



3. Stray-light correction 

The Can H channel of POLIS has two major contributions of 
stray-light (BEl 1): a spectrally undispersed ("parasitic") part,yS, 
caused by scattering inside of the spectrograph itself, and spec- 
trally resolved stray-light, a, caused by optics inside and up- 
front of POLIS. For the current investigation, we used a sim- 
plified version of the stray-light correction described in BEll. 
The observed spectra after subtraction of the dark current and 
flat fielding, /j-aw, were first corrected for the prefilter transmis- 
sion curve and normalized on average to the Liege spectral atlas 
(Del bouille et al. 1973) in the line wing near 396.5 nm (vertical 
black line in Fig.[T]). The normalized spectra /raw,norm. were then 
corrected for stray-light by 



= (/re 



- 0.025) • 0.85 . 



(1) 
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X [nm] 

Fig. 1. Intensity normalization and stray-light correction of the 
average Can H profile. Thick black: LIEGE atlas profile. Thin 
red: normalized average observed profile /raw,norm.. Thick or- 
ange dash-dotted: reference profile of the modified HSRA 
model. Thin blue: stray-light corrected average observed profile 
^orr.,norm. • The black Vertical line denotes the wavelength of in- 
tensity normalization. 

This corresponds to assuming a parasitic stray-light level of ft = 
2.5 % of Ic and a spectrally resolved stray-light level of a = 
15 %. The former value is identical to that derived in BEl 1 (5 % 
of the intensity at 396.4 nm ~ 2.5 % of 7^), whereas a is slightly 
lower. 

The values of a and JS used here were determined to match 
the average observed profile and a reference profile related 
to the inversion of the spectra. The reference profile (dash- 
dotted orange line in Fig. [T]) corresponds to the LTE syn- 
thesis of the Can H spectrum for a model atmosphere with- 
out a chromospheric temperature rise constructed from a com- 
bination of the Harvard Sm ithsonian Reference Atmosphere 
(HSRA, Gingerich et al. 1971) and theHolweger-Mueller model 
(HOLMUL, Holweger & MuUe] [19741). The exact shape of 
the reference profile in the line-core, i.e., from 396.83 nm to 
396.88 nm, where the emission peaks and the central absorp- 
tion core are located, is uncritical because only the line wings 
are matched by the stray-light correction, not the very line core. 
Therefore the choice of which temperature stratification to use 
for creating the reference profile (HOLMU L, HSRA, VAL, FAL; 
IVernazza et al.lll98"H iFontenla et al.l 1200 6*) is not important be- 
cause all these theoretical atmosphere models only diff'er signif- 
icantly in chromospheric layers, and hence in the shape of the 
line core, whereas the line wing is similar for spectra synthe- 
sized from any of the models. 

4. LTE inversion 

Initial attempts to fit the Ca spectra with the SIR code (BE08, 
Appendix B) showed that the code tends to ignore the line-core 
region in the fit. This is caused by the equal weight used for 
every wavelength point, which is perfectly suitable for photo- 
spheric lines, but not for chromospheric ones. The line core and 
the neighboring wavelength points (±0.05 nm = ± 25 pixels of 
1.92 pm in the POLIS spectra) contain all the information about 
the "upper" part of the atmosphere (log r < -2), whereas the re- 
maining about 270 wavelength points contain (redundant) infor- 
mation about the lower atmosphere. The intensity near the line 
core is the lowest of the full spectrum, and therefore any mis- 
matches between synthetic and observed spectra stay small and 
have no strong eff'ect on the least- square value to be minimized 
in the fit. 



A second problem is that the temperature stratifications re- 
quired to reproduce the observed spectra may have very com- 
plicated shapes such that an iterative modification of a given 
simple initial model may never result in the necessary com- 
plex shape, regardless of the number of nodes used for modi- 
fying the temperature stratification during the fit. To overcome 
this limitation in the inversion of the Can H spectra, we used 
as first step a comparison of the observed spectra to a pre- 
calculated archive of Ca ii H spectra, giving an increased weight 
to the line core at this step. Similar archives have been used 
in inversion appr oaches of other chromosphe r ic lines before. 



mainly for Ha ftM olowny-Horas et al. 1999|;_ 



mamiy 
l200TI: [5 



Tziotziou et al.l 



Schmieder et al...2003:.BerHcki et al...2005l) . 



4.1. Creation of an archive of LTE spectra 

The base of the archive is the modified version of the HSRA that 
emulates an atmosphere in radiative equilibrium. The archive 
profiles were then generated by adding temperature perturba- 
tions of varying shapes and amplitudes to the modified HSRA 
model and synthesizing the resulting spectra with the SIR code. 
SIR uses complete frequency redistribution and assumes LTE 
conditions. We only modified the temperature and kept all other 
parameters as given in the modified HSRA model. 

One example of the temperature perturbations applied in cre- 
ating the archive was of Gaussian shape, with diff'erent central 
positions P in log r, widths cr, and amplitudes Tampi (at log t=0). 
The ranges for the three parameters are given in Table [T] The 
temperature amplitude for the Gaussian was scaled up with its 
location in optical depth by the square-root of the electron pres- 
sure, normalized to its value at log t=0 (cf. the top panel of 
Fig. [21). The Gaussian perturbation was then moved across all 
75 grid points in optical depth from log T=-hl.4 to -6. An ex- 
ample of the resulting temperature stratifications for one run of 
the Gaussian through optical depth is given in the middle panel 
of Fig. [21 The lower panel shows the resulting synthetic spec- 
tra. This setup was chosen because spectra from a single run in 
optical depth roughly reproduce the temporal evolution of ob- 
served spectra during the occurrence of a bright grain (see the 
right panel of Fig. [4l later on). Note that because of the scaling 
of the temperature amplitude with optical depth, an initial am- 
plitude To at log T=0 of, e.g., lOOK converts to a perturbation 
with an amplitude of about 4000 K when the Gaussian is located 
at log T ~ -4. The archive thus automatically covers a large 
range of possible temperature values. We additionally globally 
added (subtracted) up to ±300 K to (from) all depth points (Tdi/f, 
first column of Table [B when running the Gaussian perturbation 
through optical depth. Additional synthetic spectra were created 
by adding (subtracting) straight lines of variable slope to (from) 
the modified HSRA model instead of adding a Gaussian pertur- 
bation. 

Line-of-sight (LOS) velocities were not included in the gen- 
eration of the archive profiles because this would have increased 
the size of the archive beyond usefulness. Velocities are instead 
taken into account in an empirical way in a later step. The spec- 
tral resolution of the Ca ii H spectra from POLIS corresponds to 
a velocity dispersion of about 1.5 kms"^ per pixel. Together with 
the intrinsic (thermal) width of spectral features, the weak ve- 
locity response allowed us to still use the pre-calculated archive 
spectra as first step to fix the temperature stratification because 
the Doppler shifts (in pixel) stay comparably small. The macro- 
turbulent velocity in the atmosphere was set to a fixed value of 
1 kms"^ in the generation of the archive spectra, which yielded 
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Fig. 2. Top: scaling law for the temperature amplitude with opti- 
cal depth. Middle: example of temperature stratifications with 
the Gaussian perturbation centred subsequently at decreasing 
optical depth. Bottom: the 75 synthesized spectra corresponding 
to the temperature stratifications above for the Gaussian pertur- 
bation centred subsequently at the location P from -\-lA (bottom- 
most spectrum) to -6 (topmost spectrum) in log r. 

Table 1. Parameter ranges used in generating the archive. The 
values in parentheses give the step width. 

TdifflK] Ta^pi\r=i [K] o- [logT] P [logr] 

-300 - +300 (50) 10- 100(10) 0.3 -3.1 (0.1) -6 -+1.4(0.1) 

line widths for the photospheric line blends that match the ob- 
served line widths. 

With the Gaussian perturbations defined by the free param- 
eters (Tdiff, P, cr, Tampi) and the other additional modifications 
of the temperature stratification in total more than 300.000 Ca ii 
H spectra were synthesized with SIR. All line blends of Can 
H with known transition parameters were taken into account 
(cf. Table 1 in BE09). We also always synthesized from the at- 
mospheric models the corresponding spectra at 630 nm because 
POLIS provides both wavelength ranges simultaneously. The 
creation of the archive by a single (background) process takes 
about one week. 

4.2. Best-fit spectrum from LTE arcliive 

To find the archive profile A(A, i) (i counts the number inside 
the archive) that matches best a given observed profile 1(A), 
the least-square deviation between the observed profile and all 
archive profiles was calculated by 

j 



The weights w(Aj) were set to the squared inverse of an average 
QS profile, w = 1 //^^, with a subsequently doubled weight for 
the region around the line core itself from about 396.81 nm to 
396.89 nm. This strongly enhances the influence of the line-core 
intensities relative to the wing (ratio of about 4:1). All (strong) 
line blends were masked out by setting the weight to zero at these 
wavelengths. The 630 nm spectra were not included in the;^^ for 
finding the best-matching archive profile because the fit to the 
Ca line core is the main driver for this step. 

4.3. Iterative improvement of temperature stratification 

The best match of observed and archive profile is determined 
with a strong weight for the very line core. In the line wing, 
the observed and best-fit archive profile usually deviate because 
the archive only contains a limited number of all of the possible 
combinations of relative wing and core intensities. To improve 
the fit in the line wing, we used an iterative modification of the 
temperature stratification based on the intensity diff'erences be- 
tween observed and synthetic profiles at some wavelengths. The 
temperature values at the corresponding layers in log r that con- 
tribute to the intensity at these wavelengths are then modified 
according to the modulus and sign of the intensity diff'erences. 

We selected a series of wavelength windows without photo- 
spheric blends and the blue part of the line core up to the line 
centre (cf. the black crosses in the top row of Fig. [3]). For these 
wavelengths A, we calculated the corresponding layers in log r 
using the intensity response function of Ca ii H in the modified 
HSRA model (cf. iRezaei etal.ll2QQ8l BE09). The intensity re- 
sponse function matches well formation heights derived from 
phase diff'erences of propagating waves (BE08, BE09). We then 
implemented the relation between A and corresponding r into the 
code (provided by an external file). It turned out that the hard- 
wired response function works fine for all spectra on disc centre, 
but for the AR observations at a heliocentric angle of 50°, the re- 
lation between A and r was already slightly off' leading to a mis- 
fit in the outermost wing/lowermost atmosphere layers. It seems 
therefore to be recommended to at least re-calculate the average 
A-T relation for an application to off'-centre observations. 

During the iterative improvement of the complete temper- 
ature stratification also the simultaneous 630 nm spectra were 
now included by adding one more point to the wavelength win- 
dows used. This point was located in the continuum of the 
630 nm channel, corresponding to logr ~ 0.2 in the intensity 
response. The lowermost forming wavelengths in the Ca wing 

in the usual setup of POLIS form at about logr 0.2 (cf. the 

black crosses in the middle row of Fig. [3]) The use of the 630 nm 
continuum intensity was, however, implemented only optionally 
and can be skipped when no 630 nm data are available. 

With the relation between A and r given, one can then mod- 
ify the temperature stratification according to the intensity dif- 
ferences between observed and best-fit spectrum, A/(/l), at the 
wavelengths considered. The intensity diff'erence is first con- 
verted to a temperature diff'erence using the intensity response 
function. We then fit a fourth-order polynomial to the required 
temperature changes in log r to obtain a smooth curve for all op- 
tical depth values. The slope of the curve is extrapolated towards 
the layers with log r > 0.2 where no measurements are available. 
The change of the temperature stratification is then added to the 
previous stratification, and the new spectrum is synthesized. The 
middle panels of Fig. [3] show some examples of the variations 
during the iteration of the temperature stratification. The method 
converges usually on two, three iterations. The code uses no spe- 
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cial convergence criterion, only the number of iterations to be 
done has to be provided by the user. 

In the iteration, no special weight is given to the line core 
because this step is mainly intended to improve the fit in the line 
wing. As a result, the;^^ value using the weighting with an en- 
hanced contribution of the line core (cf . the previous section) can 
actually worsen. Even if this happens, the iterated temperature 
stratification is returned because otherwise it is not possible to 
obtain a good fit to the line wing. To identify such cases, i.e., 
a degradation of the fit to the Ca line core in the iteration pro- 
cess, the inversion routine returns an additional binary value of 
or 1 that indicates a decrease or increase of the weighted 
(cf. Sect. [O later on). 

The LTE condition in the spectral synthesis of the profiles 
requires the temperature stratifications to have a negative slope 
in the uppermost layers for creating double reversals in the line 
core (cf. the middle right panel of Fig. O. This usually aff'ects 
the last 3-5 grid points of the optical depth scale (logr = -5.5 
to -6) at the upper end of the optical depth scale. The profiles 
created by the addition of the Gaussian perturbation to the mod- 
ified HSRA automatically have such a negative slope whenever 
the centre of the Gaussian is not located exactly at log r = -6. 

4.4. Inclusion of LOS velocities 

In the iterative improvement of the temperature stratification, the 
LOS velocities are finally included in an empirical way. We use 
three of the line blends of Can H (Cri at 396.37 nm, Fei at 
396.61 nm, Fei at 396.93 nm) and the Ca line core itself. The 
line blends in the wing correspond to three diff'erent formation 
heights (cf. BE09). We use the observed LOS line-core velocities 
of the four lines to construct a velocity stratification by attribut- 
ing the observed velocities to the corresponding log r layers of 
-2.4, -3.4, -3.9, and -6, and then linearly interpolate between 
the velocity values. To prevent step functions, the resulting ve- 
locity curve with optical depth is smoothed. The bottom panels 
of Fig. [3] show the resulting velocity stratification for the ob- 
served spectra in the top panels. The velocity is not modified to 



Fig. 3. Iterative improvement of the 
best-fit archive spectrum for two exam- 
ple spectra without (left column) and 
with emission in the line core (right 
column). Top row. observed spectrum 
(black line), best-fit archive spectrum 
(blue), and final best-fit spectrum (red). 
The black crosses denote the wave- 
lengths and intensities used for the im- 
provement of the fit. 2nd row. temper- 
ature stratifications (blue: archive best- 
fit, red: final best-fit, black: modified 
HSRA model). The black crosses de- 
note the optical depths corresponding 
to the wavelengths marked above. The 
diamond at log r =0.2 denotes the 
630 nm continuum intensity contribu- 
tion. Bottom row: stratification of the 
LOS velocity. 



obtain a better fit, it serves only to yield a better match of syn- 
thetic and observed spectra using the observed velocity values. 
For the line blends, the velocity is well determined, whereas for 
the Ca line core the location of the line minimum in some cases 
might not reflect a true velocity at all. 

For the inversions of the AR maps, we did not include LOS 
velocities because the line core of Ca ii H has there often a com- 
plex shape with a single (or several) intensity reversals that does 
not easily allow one to ascribe any velocity to it. Because of the 
weak response of the spectra to LOS velocities, this should have 
had a negligible efl'ect on the retrieved temperatures. 

It takes about 8 seconds to obtain a fit to a single observed 
profile on a common desktop computer, and hence about three 
days for all observed spectra in one typical observation of 30.000 
profiles. 



5. Results 

5.1. Example spectra in quiet Sun 

For a visual inspection of the quality of the fits, we show a se- 
ries of spectra in the following Figs. IH to [HI The left panel of 
Fig. |4] shows the observed and best-fit Ca spectra (first and sec- 
ond column) for three locations along the slit in the time- series 
in QS. The third column shows the corresponding temperature 
stratifications. The fourth column shows the simultaneous and 
co-spatial Stokes V spectra at 630 nm to demonstrate the pres- 
ence (or absence) of photospheric magnetic fields. The upper 
row correspon ds to a location in or near the photospheric net- 
work (see also iBeck et al.ll2008ab . An animation of the spectra 
along the slit and the corresponding inversion results for all time 
steps is available in the online material (see Appendix lAb. 

Comparing observed and best-fit spectra in Figs.lHtoO we 
find that the inversion approach is able to deal well with the spec- 
tral range from the wing up to about 396.8 nm, fully reproducing 
the intensity patterns and the Doppler shifts of the line blends. In 
the Ca core, the red emission peak is in general always less well 
reproduced than the blue one. This is caused by the inability to 
reproduce the line-core region without including NLTE eff'ects 
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Fig. 4. Comparison of observed spectra and inversion results in the time-series. Left panel, top to bottom: spectra from three different 
locations along the slit versus (vs.) time. The upper row corresponds to a network location. First column: observed spectra. Second 
column: best-fit spectra. The third column shows the temperature stratifications vs. logr, the fourth column the simultaneous Stokes 
V spectra of the 630 nm channel. The white rectangle in the second row marks the section shown magnified at the right-hand side. 
A train of three consecutive brightenings passes from the line wing towards the line core between t = 300 sees and t = 800 sees. 
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Fig. 5. Same as Fig.lHfor the large-area scan. Note that this time the cuts are along a spatial axis and do not represent a temporal 
evolution. The white rectangle in the top row marks the section shown magnified at the right-hand side. 
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Fig. 6. Example spectra during the passage of a 
(shock) wave. Observations are given by black 
pluses and the best-fit spectra by red lines. The 
temporal cadence between subsequent spectra 
is 21 sec. Time increases from left to right in 
each row and from bottom to top between rows. 
The coloured bars in the upper left panel de- 
noted the wavelength bands used later on. 
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Fig. 7. Same as Fig.[5]for three cuts through the 
active region. No LOS velocities were consid- 
ered in this fit. The temperature is displayed on 
a logarithmic scale in a range from 3500 K to 
10000 K this time. 
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and mor e complex velocity fields than used in the inversion (cf . 
lCMsson& Stein 1997). 

The core of the spectral line of Fei at 396.93 nm reverts to 
emission in the best-fit spectra on several locations, contrary to 
the observations (e.g., at ^ = 400, 600, 750 sees in the right panel 
of Fig. (4]). The LTE assumption thus already fails for the core of 
this line in some cases. The temporal evolution, e.g., wave trains 
propagating from the line wing towards the core, and all vari- 
ations in the line wing are well captured by the inversion. The 
behaviour of the observed spectra during the passage of three 
consecutive wave fronts {right panel of Fig. |4]) is well matched 
except for the asymmetry of the red and blue emission peak. The 
corresponding temperature stratifications show both the upwards 
propagation of the wave fronts and the increase of their ampli- 
tude in the upper atmosphere. 

The diflTerences between magnetic and field-free locations 
can be seen better in Fig. [5] that shows spatial cuts through the 
large-area map. Whereas the (temporary) enhancements typi- 
cal for propagating waves often only appear in the upper lay- 
ers (log T < -4) with a reduction of temperature between about 
log T = -2 to -4 (cf . the top panel at the right-hand side of 
Fig. 131), the temperature on most locations with significant po- 
larisation signal, and hence magnetic fields is usually increased 
over the complete optical depth range {third column of the left 
panel of Fig. O. The magnification at the right-hand side of 
Fig. [5] shows that in some cases regions with reduced inten- 
sity in the very line core can extend spatially connected over 
5-10'', reflected by the corresponding reduction of temperature 
diiy = 20 - 30'', in contrast to the magnetic location around 
y = 9". 

The accuracy in the reproduction of the observed profiles in 
the QS is visualized in more detail in Fig. [6l that shows the evo- 
lution of individual spectra during about three minutes in the 
time-series. In that time span, one wave front passes through the 
Ca line formation range, culminating in strong emission of the 
H2V peak when reaching the line core. The intensity near the line 
core rises gradually (profiles Nos. 23 and 24), yielding a profile 
with an elevated line wing close to the line core and a narrow 
absorption core without emission peaks in profile No. 25 {left- 
most middle panel of Fig. [6]) that gives way to one with a strong 
H2V emission peak 20-40 sees later at profile No. 27 {rightmost 
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Fig. 8. Individual spectra from the active region map. Clockwise, 
starting left top: quiet Sun, bright point in QS, umbra, penumbra. 
Black: observed spectrum. Red: best-fit spectrum. 



middle panel of Fig. [6]). The emission in the line core subsides 
again after the passage of the shock front (profiles Nos. 28 to 
30), whereas at the end in the line wing already the next wave 
front is indicated (compare profiles Nos. 24 and 30). The indi- 
vidual profiles clearly demonstrate that the H2R emission peak is 
not reproduced by the fit, but the amplitude of H2V is matched 
well, e.g., in profile Nos. 26 to 29. 

5.2. Example spectra in active region off disc centre 

Figure |7] shows observed and best-fit spectra for three spatial 
cuts through one of the three active region maps in the same 
layout as in the left panel of Fig. |4l The best-fit spectra diflTer 
more from the observed spectra than for the QS data on disc 
centre. No LOS velocities were included in the fit because for 
the line core of Ca 11 H it is difl&cult to determine any suitable 
value directly from the spectra. The profiles in the umbra show 
only a single emission peak (where a fit of a Gaussian could 
be used), those in the penumbra usually double reversals (where 
the location of the central absorption could be used), but also 
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Fig. 9. Bottom panel: spatio- 
temporal maps of wavelengths 
bands in the observed spectra {odd 
rows) and in the best-fit spectra 
{even rows) of the time series. The 
lower two rows show from left 
to right the continuum intensity 
at 630.4 nm, OW, MWl, MW2, 
MW3, and IW. The upper two rows 
show IWl, H2V, core intensity, 
H2R, RWl, and RW2. Top panel: 
corresponding temperature maps 
at several optical depths from -0.5 
to -5.5 in logr (denoted at the top 
in each subpanel). The minimum 
and maximum temperatures are 
denoted at the left and right end 
of the grey bar, respectively. The 
three values on each side in the 
lower row correspond to the layers 
of logr = -0.5/ - 1.5/ - 2.5 (from 
left to right), those in the upper row 
to logr = -3.5/ -4.5/ -5.5. 



mixtures of these two types or profiles with a plateau of emission 
appear. A determination of LOS velocities to be attributed to the 
Ca line core would thus require a (predetermined) choice of the 
method to derive the velocity because it is impossible to use a 
unique method over the complete FOV. However, because of the 
small impact of LOS velocities up to the sound speed on the 
spectra, neglecting the post-facto inclusion of velocities should 
not have influenced the retrieved temperature stratifications. 

Nevertheless, a clear deviation between best-fit and observed 
profiles is seen in the very line core of Ca ii H. All best-fit pro- 
files outside the umbra are characterized by a central absorp- 
tion core that is neither seen in the observed AR spectra nor the 
QS data. This fact is better visible in the individual spectra of 
Fig. [8] that correspond to a position in the umbra, penumbra, a 
bright point (BP) in a plage area, and a "QS" location in the AR. 
Except for the umbral profile, all best-fit profiles have a dark 
central absorption core. This absorption core presumably results 
from the requirement to have a decrease of the temperature in 
the uppermost layers to obtain " reasonable" s pectra (imposed 
by the LTE assumption, cf. also iRezaei et al.i r2008). It shows 
up much more prominent for the AR LTE spectra because there 



a relatively strong chromospheric temperature rise is required 
to produce the strong emission core at first, on whose top the 
(small) absorption core then is superposed. In the QS spectra, 
the strong emission core is absent. Another diff'erence between 
Figs. [5] and [7] is the temperature range in the AR required for 
displaying the temperature stratifications in the third column of 
Fig. [71 The variations of temperature are so much larger than in 
QS that the display range had to be significantly increased, even 
while the display now is on a logarithmic scale. 

The creation of the archive spectra was done in a semi- 
automatic way without explicitly taking into account diff'erences 
between spectra in QS and AR. Figure [8] demonstrates that the 
archive by chance also covered the range of profile shapes that 
are observed in ARs. 



5.3. 2D maps in quiet Sun 

The lower panel of Fig. [9] shows 2D spatio-temporal maps of 
the intensity in selected wavelength bands as final cross-check 
of the performance across the comp lete FOV and the full spec- 
tral range (compare to, e.g.. iCauzzi et al.ll2008bl their Figs. 2 
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Fig. 10. Bottom panel. 2D maps of wavelengths bands in the observed spectra {odd rows) and in the best-fit spectra {even rows) of 
the large-area scan on disc centre. The lower two rows show from left to right the continuum intensity at 630.4 nm, OW, MWl, 
MW2, MW3, and IW. The upper two rows show IWl, H2V, the line-core intensity, H2R, RWl, and RW2. Top panel, corresponding 
temperature maps at several optical depths from -0.5 to -5.5 in log r (denoted at the top in each subpanel). See Fig.[9]for additional 
information on the grey bar. 



and 3). The wavele ngth bands are described in more detail in 
iRezaei et"aD (l2007 h and BEOS, and are marked in the upper left 
panel of Fig.O They pass from the outer line wing (OW) in the 
blue through the line core towards the red line wing (RW), avoid- 
ing the line blends. The only diff'erences between observed and 
inverted spectra that can be discerned by eye are in the map of 
the continuum intensity at 630.4 nm and in the H2R map. For the 
former, the dynamic range of the observed continuum intensity 
at 630.4 nm is not fully matched and the spatial pattern in the 
inversion spectra is more similar to that in the OW than to the 
observed continuum intensity itself. This presumably is related 
to the small influence of the continuum intensity in the itera- 
tive improvement of the fit because it provides only one wave- 
length, and hence one optical depth point to the modification of 
the temperature stratification (Sect. 14.31) . For H2R, the observa- 



tions show a higher intensity than the fit for the network regions 
at 3; - 10 - 20'' and around 50''. 

The upper panel of Fig. [9] shows temperature maps at sev- 
eral optical depth values for the inversion of the time- series, for 
a visualization of the relation between observed spectra and the 
corresponding temperatures. The LTE assumption maps the in- 
tensity at some wavelength to the temperature at some optical 
depth, but because the radiative transfer through the complete 
depth stratification was considered, the approach exceeds the di- 
rect conversion from intensity at one wavelength to temperatures 
through the Planck function used in, e.g., Cauzzi et'aP (l2009h or 
BE 12. During the iteration of the best-fit stratification from the 
LTE archive the code also always modifies the temperature over 
some range in log r such that the profile after the integration of 
the radiative transfer equation matches the observed spectrum. 
This also breaks up the one-to-one correlation between a given 
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Fig. 11. Temperature maps for the active region scans. Note that the layers of log r slightly differ (0 to -5 in log r). The temperature 
values of the grey bar refer to the middle row as reference, the range in the other two scans is of comparable order. 



wavelength and a single layer in log r. Comparing the intensity 
and temperature maps in Fig. |9j one can identify some features 
that appear modified in the spectra and the temperature strati- 
fication. One case is a large-scale c-shaped darkening around 
it = 40min, y = 25'') that appears only prominently in the 
line-core map of the observations and weakly in H2V and H2R, 
but can be seen in the temperature maps for all layers above 
logr = -4.5. 

Figure [To] shows the intensity maps in the wavelength bands 
in the large-area scan in QS on disc centre. For the compari- 
son between observed and best-fit spectra in Fig. [TOl the match 
seems to be even better than for the spectra of the time- series 
on disc centre shown in Fig. |9l The increasing spatial extent of 
the emission near the locations of photospheric network fields 
in wavelength bands near the line core is clearly seen in both 
observed and best-fit spectra. Wavelength bands that form lower 
are also well matched. The largest deviations can be seen in the 
map of the H2R emission peak (fourth column in the upper two 
rows), the RWl map (fifth column in the upper two rows), and 
the continuum intensity at 630.4 nm. The latter two maps show 
some intensity oflfset and contrast diff'erence between observa- 
tions and best-fit spectra. 

The temperature maps in the upper panel of Fig. [TOl show the 
extension of the emission around the locations of concentrated 
photospheric magnetic flux for logr smaller than -4.5. Isolated 
temperature enhancements (e.g., at x,y ~ 65'', 50'') seem to be 
more pronounced and localized than the corresponding brighten- 
ings in intensity in H2V or the line-core map. The magnetic loca- 
tions are not discernible in the temperature map at log r = -1.5. 
The lowest layer displayed at log r = -0.5 already clearly forms 
above the continuum layers and corresponds roughly to the in- 
tensity map of MW2. 

5.4. 2D maps in active region 

In the temperature maps of the AR (Fig. [TT]) several peculiari- 
ties stand out. The map shown in the bottom row exhibits two 



dark filaments for logr < -4 at x ~ 70" and y ~ 15" and 
30", respectively. The lower, larger one has a co- spatial filament 
in He I spectra at 108 3 nm, and was identified with a transient 
siphon flow by Bethge'etal] (l2012h . Several similar dark fila- 
ments that originate or terminate in individual pores can be iden- 
tified in the scan shown in the middle row. An investigation of 
the related Doppler shifts will allow one to determine if they are 
also related to si milar siphon flows (cf. ^ Uitenbroek et al.lE006l: 
iBeck et aLl lMol). The AR map in the top row shows two inter- 
esting emission features, or equivalently locations of strong tem- 
perature enhancements. The first is an isolated bright point (BP) 
at (x, y) ~ (45", 30"). This BP has about the highest temperature 
of the full FOV at logr = -3, and still stands out in the map 
at logr = -4. The Ca spectra at this locati o n are similar to the 
BP spectra shown in Fig. 10 of iBeck et al.l (120051) . The second 
large-scale region of temperature enhancements is located just 
between the large and small sunspot in the top row of Fig. [TT] 
at (x,y) ~ (40", 10" - 20"). The increase of emission between 
the two sunspots at log r = -3 to -5 could indicate the pres- 
ence of a current sheet that s eparates their field line s, similar to 
the current sheet discussed in Tritschl er et al. ' (I2008h for Ca II IR 
observations (see also lGoodman & Judge 20121) . 



6. Discussion 

6.1. Quality in reproduction of observed spectra 

The comparison of the 2D intensity maps of observed and best- 
fit spectra showed some diff'erences between the two. To quantify 
the resulting mismatch, we made scatter plots of the observed 
and best-fit intensities in the wavelength bands for the large-area 
scan on disc centre (Fig.[T2|. The graphs confirm that the inten- 
sity in the continuum at 630.4 nm and in the H2R emission peak 
deviate between observations and best-fit spectra: for the for- 
mer, the slope diff'ers from unity, while for the latter the scatter 
is significantly larger than in, e.g., H2V. For all other wavelength 
bands the 3-cr interval of scatter (red lines) covers the line of 
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Fig. 12. Scatter plots of observed and best-fit 
intensities {black dots) of the large-area scan on 
disc centre. The blue dots denote binned values 
whose 3-cr variation is denoted by the red lines. 
The solid black line denotes a one-to-one cor- 
relation. Bottom right panel: standard deviation 
of the diff'erence between observed and best-fit 
intensity as function of wavelength. 



one-to-one correlation. The bottom right panel of Fig. [T2l shows 
the standard deviation of the difference between the observed 
and best-fit intensities as a function of wavelength. Besides from 
the location of the line blends, the rms deviation is about 0.5 % 
of Ic in the line wing, and increases to about 1 % of for wave- 
lengths in the line core from 396.82 to 396.88 nm. The rms at the 
location of the H2R emission peak is slightly enhanced relative 
to rest of the line core. 

The spatial variation in the reproduction of the observed 
spectra is displayed in Fig. [131 In the map of with equal 
weighting for all wavelengths, i.e., w{X) = 1, that measures 
the quality of the fit to the full line profile {top right panel), 
the centres of the network regions stand out, but only slightly. 
The cellular granular- sized pattern in the;^^-map corresponds to 
spatial structures seen in wavelengths forming at and below the 
MW and indicates misfits in the lower atmosphere. Calculating 
the with the enhanced weighting for the line core as used 
in the retrieval of the best-fit archive profile {bottom left panel) 
yields a similar spatial variation of the fit quality, with the net- 
work locations now standing out more prominently and a more 
homogeneous quality throughout the inter-network regions. The 
marker that shows how the x^ with the enhanced weighting for 
the line core was modified by the iterative improvement of the fit 
is displayed at the lower right. It turns out that on about 60 % of 
the FOV, the iterative modification of the full temperature strati- 
fication has worsened the initial fit to the very line core obtained 
using only the best-fit archive profile. Only in the inter-network 
regions with usually less complex profile shapes, the iteration 
has improved not only the fit to the line wing, but also to the 
line core even with the strongly enhanced contribution of the 
line core. For the AR data, in nearly all cases (> 90 %) the fit 
to the line core was worsened in the iterative modification. We 
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Fig. 13. Spatial variation of the x^ in the large-area QS scan. 
Top row, left: line-core intensity as indicator of the solar spa- 
tial structuring. Top right, x^ with equal weight for all wave- 
lengths. Bottom left: x^ with enhanced weighting for the line 
core. Bottom right, decrease/increase of x^ {white/black) by 
the iterative modification of the fit when considering enhanced 
weighting in the line core. 



note that this is a slight worsening relative to the best-fit archive 
profile, with still an acceptable fit to the line shape even for most 
of the AR spectra (e.g.. Fig. (8]). 
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Fig. 14. Relative occurrence of a given archive profile in the in- 
version. Thick black: QS map. Thin red: AR map. 



6.2. Statistics on the arcliive usage 

The size of the pre-calculated spectral archive partly limits the 
usefulness of the two-step inversion approach described here. At 
the spectral sampling of the POLIS data of 1.92 pm per pixel 
and the total wavelength range covered, the full archive of about 
300.000 spectra could be easily kept in the memory of a sim- 
ple standard desktop machine because its size was only about 
350 MB. Figure [T4l shows the statistics on how often which 
archive profile was actually used in the determination of the best- 
fit archive profile. The most frequently used archive profiles ap- 
pear with about 3 - 4 % relative frequency. Some areas in the 
archive were only relevant in either the QS or the AR data, e.g., 
the archive profiles with numbers above 200.000 only appear for 
the AR spectra. These profiles correspond to a large positive off'- 
set Tdiff in temperature in their synthesis. This global off'set in 
temperature between the QS and AR data is caused by both the 
change in the heliocentric angle, where in the AR data the LOS 
reaches a given optical depth at a diff'erent geometrical height in 
comparison to the QS data, and the presence of magnetic fields 
throughout nearly all of the FOV in t he AR that is known to lead 
to i ncreased emission in Ca spectra (iRezaei et al .1120071: iRezaeil 
12008). 

It turns out that the actual usage of the archive is much more 
sparse than the plot of Fig. [14] suggests. In total, only about 5700 
diff'erent archive profiles (~ 2 % of the full archive) were picked 
in the inversion of either QS or AR data (-80.000 observed spec- 
tra). The acceptable quality of the fits in both QS and AR implies 
that an exchange of the LTE archive for about 10.000 profiles 
calculated in NLTE conditions could already suffice to provide a 
reasonable NLTE fit to observed chromospheric Ca ii H profiles, 
and hence as well to similar spectra in other chromospheric lines. 
The full dynamical range of the temperature variations seems to 
be more critical in the creation of the spectral archive than the 
density of the sampling inside the parameter space of tempera- 
ture. 



6.3. Outiook 

Our findings can be extended to what will be needed for the fu- 
ture analysis of spectra from, e.g., the planned Solar-C mission, 
regardl ess of a decision for the Ca ii H and K or the Mg ii h and 
k lines (iKatsukawa et al.ll201 ll:lBelluzzi & Truiillo Buenoll2012h 
because of their similar behaviour, or data from the p lanned 
BLue Imaging Solar Spectrometer (BLISS, Puschman n et al.l 
2012bl in press) fo r the new 1.5-m GREGOR telescope (e.g., 
Schmidt et al.ll2012l in press). The BLISS data will be accompa- 



nied by spectropolarimetric observations at, e.g., 630 nm with 
the G REGOR Fabry-Perot Interferometer (Puschmann et al.l 
l2012al in press), which will provide the possibility to simulta- 
neously measure photospheric magnetic fields as in the case of 
POLIS. 

For deriving the thermodynamic properties of the solar chro- 
mosphere from observed spectra, direct analysis tools will be 
few and very limited. The best approach seems to consist in a 
"controlled" NLTE forward modeling, i.e., a creation of NLTE 
spectra for what can be considered realistic solar atmosphere 
models. These could be derived fro m both analytica l 2D so- 
lutions of solar surface structures (Uitenbroek 2011) or syn- 
thetic spectra from numerical simulations (L eena arts et al. 2009j, 
120101: IWedemever-Bohm& Carlssonll201lh . The lateral radia- 
tive transport is more important for chromospheric than for pho- 
tospheric lines. A creation of an archive from 3D numerical sim- 
ulations might thus be anyway r equired for a consistent spectral 
synthesis (cf. Uitenbroek & Criscuolil 1201 ll : ISheminoval |20 1 21: 
iLeenaarts et al.ll2012h . Such an archive then also allows one to 
improve the inversion procedure beyond the common single- 
pixel approach in which every pixel is inverted individually ne- 
glecting its surroundings. Especially for chromospheric lines, a 
simultaneous fit to a small-scale area of, e.g., about x (3 x 3 
pixels at the POLIS spatial sampling) can be crucial to obtain 
consistent results. An NLTE archive from 3D numerical simula- 
tion can directly provide suited spectra on a 2D area including a 
central pixel and its surroundings. There will be a need for addi- 
tional modification of initial profiles and the related thermody- 
namical variables to match the observed spectra, similar to the 
two-stage approach taken here, because the combined parame- 
ter space of the uncorrected temperature and velocity variations 
cannot be covered by any reasonably sized archive. 

The creation of a base of NLTE spectra seems to be a promis- 
ing approach for both the Ca and Mg lines, with the advantage 
that such spectra can already be creat ed and tes ted against si mu- 
lation s and/or theoretical models (cf. de la Cruz Rodriguez et al.l 
[2QT2h well in advance of the real data to be obtained in a few 
years. 

The results of the retrieved temperature stratification in QS 
and AR will be discussed in the subsequent paper of this series. 

7. Conclusion 

A pre-calculated archive of about 300.000 spectra is found to 
be more than sufficient to fit the line cores of observed Ca ii H 
profiles in quiet Sun and active regions satisfactorily regardless 
of the huge range of possible shapes, whereas the line wing has 
to be iteratively modified, because line core and line wing vary 
independently of each other. The advantage of the two-step in- 
version approach is that an exchange of the archive for a suited 
data base of spectra calculated in NLTE suffices to provide a 
first-order NLTE inversion method. 

We have to leave it up to the reade r to decide if we met the 
requirement given in the quotation from'Cram' (1911) at the very 
beginning to provide a "careful study". By using the LTE as- 
sumption, we are, however, perfectly sure that we fulfill his la- 
beling of using a "crude" method. 
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Appendix A: Animation of inversion results 

Figure IAJI shows a still from an animation that displays the in- 
version results and some line characteristics from both the Ca 
and the 630 nm spectra for the time-series of quiet Sun. There 
are three 2D maps: the H-index (integration of the line-core in- 
tensity of Ca) in the first panel, the polarity of the Stokes V signal 
(±1 for magnetic fields along/anti-parallel to the LOS and zero 
for locations without significant polarization signal) in the sixth 
panel, and the LOS velocity of the Fei line at 630.15 nm in the 
seventh panel. For all these 2D maps, eight repetitions of scan- 
ning the same area of x 61'' are shown, i.e., the temporal evo- 
lution inside the area over the course of about 170 seconds is dis- 
played. The spectra shown in the rest of the panels were always 
taken from the scan step that is located at 0'' in the 2D maps. The 
animation runs over the complete duration of the time-series of 
about one hour. 

In the 2D maps of the H-index, the lateral spatial motion of 
individual bright grains can be followed. The temporal evolution 
of the temperature stratifications shows the spatial (along the slit 
and in height) and temporal extent of wave trains and tempera- 
ture perturbations. One prominent example of an upwards propa- 
gating temperature perturbation that subsequently aff'ects the full 
line profile appears between t = 378 and 483 seconds at j ~ 28'', 
i.e., at about the middle of the FOV along the slit. Lateral mo- 
tion of temperature perturbations along the slit can be seen for 
some events that happen close to locations with strong polar- 
ization signal. One such event shows up between ^ = 693 and 
756 seconds, starting at the location of the persistent polariza- 
tion signal at J ~ 40" and moving downwards along the slit with 
time, while getting increasingly weaker on the way. The Stokes 
V spectra and the polarity of the Stokes V signal show that 
most of the polarization signals, and hence photospheric mag- 
netic fields detected are rather stable throughout the time-series. 
With an integration time of less than 4 seconds per scan step, the 
polarimetric sensitivity in the time- series is naturally limited to 
strong polarization signals that correspond to the stable photo- 
spheric magnetic network. The polarity was erroneously labeled 
with "|y|" during the creation of the animation. 




Fig. A.l. Still from an animation of the inversion results for the time-series of quiet Sun. Leftmost panel. 2D maps of the H-index 
for eight repetitions of scanning the same area of 2" x 61" . Time increases left to right. The scan step of the spectra shown in 
the other panels is always located at the very left border at 0". Second and third panel: observed and best-fit Ca spectra. Fourth 
panel, temperature stratifications along the slit corresponding to the best-fit spectra. Fifth panel, co-spatial and simultaneous Stokes 
V spectra at 630 nm. Sixth panel. 2D maps of the polarity of the Stokes V signal for eight repetitions of scanning the same area. 
Seventh panel 2D maps of the LOS velocity of the Fe i line at 630. 15 nm for eight repetitions of scanning the same area. 



